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We present results of a unitary triangle fit based on the scan method. This frequentist approach 
employs Gaussian uncertainties for experimental quantities, but makes no arbitrary assumptions 
about the distribution of theoretical errors. Instead, we perform a large number of fits, scanning 
over regions of plausible theory errors for each quantity, and retain those fits meeting a specific con- 
fidence level criterion, thereby constraining the p — fj plane using the standard input measurements 
(CKM matrix elements, sin 2/3, B® s mixing, e^) as well as branching fraction and CP asymmetry 
measurements of B decays to PP,PV, VV, and a\P final states to determine a, D^K^ modes to 
determine 7, and D^tt and Dp modes to determine 2/3 + 7. We parameterize individual decay am- 
plitudes in terms of color-allowed tree, color-suppressed tree, penguin, singlet penguin, electroweak 
penguin, as well as W-exchange and IV-annihilation amplitudes. With this parameterization, we 
obtain a good fit to the measured branching fractions and CP asymmetries within the Standard 
Model ansatz, with no new physics contributions. This simultaneous fit allows us to determine the 
correlation between a and /3 as well as between 7 and /3. 

PACS numbers: 12.15.Hh, 13.25.Hw, 14.65.Fy 



I. INTRODUCTION 



The phase of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix [l[ is responsible for CP violation in the Standard 
Model (SM). The unitarity relations within the CKM ma- 
trix provide an excellent laboratory to test this predic- 
tion, the relation V ub *V ud + V cb *V cd + V tb *V t d = being 
particularly useful, since many measurements and the- 
ory inputs in the B and K systems can be combined 
for this test. Several such tests, such as those by the 
CKMfitter and XST fit groups are m common circu- 
lation. The former is a frequentist technique, the latter 
Bayesian. The scan method J4| presented herein takes a 
different approach in its treatment of the distribution of 
theoretical uncertainties. 

We first describe the scan method. We then employ 
the method to extract CKM parameters using inputs 
standardized for the book Physics of the B Factories [5| , 
comparing the results obtained with these parameters by 
CKMfitter and UT/j t , and investigating the question of 
consistency of the results with SM expectations. These 



fits utilize stand-alone determinations of the angle a and 
7 as inputs. We then extend the scan method fit to in- 
clude in the unitarity triangle parameter determination 
the measurements of B decays to PP, PV, VV, and a\P 
final states to determine a, to D^K^ final states to 
determine 7, and to D^'ir and Dp final states to deter- 
mine 2/3 + 7. I n this wa y> we are able for the first time to 
include the correlations between a and /3 and 7 and /3 in 
the extraction of the parameters of the unitarity triangle. 

II. FIT METHODOLOGY 

The scan method is a frequentist-based fitting tech- 
nique to determine the parameters of the CKM ma- 
trix. Rather than making assumptions as to the distribu- 
tion of theory errors, it accounts for the uncertainties in 
the QCD parameters f Bs , f B J f Bd , B Ba ,B Bs /B Bi , and 
Bk and the CKM parameters V ub and V cb by scanning 
over the range of theory uncertainties using fixed grid 
or Monte Carlo (MC) methods. In the baseline fit, we 
combine measurements of Amj, Am s , e^, V cb , V ub , V u d, 
V us ,sm2/3,a and 7 in the \ 2 function: 



X (p,V,Pi,tj) 



(Amj, s )-Am d|S (p,i),pi,fj) \ / {\V cb 

,ub,ud,us ,ud,us 
°"Am rf „ / I cry, 



+ 



(\eK\)-eK(p,V,Pi>tj) 



+ 



cb.ub,ud,u 

2 



(a)-a(p,n,pi) 

0~ot 



{Tn)-Tn(pi,tj) 



(1) 



where the pi are measured quantities including the parameters. To account for correlations in the different 
Wolfenstein parameters [6j A and A, and the tj are QCD 
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TABLE I. Observables used in the baseline fits. The first values of V c b an d V u b are those used for the comparison with 
CKMfitter and UT^ 4 . The second, updated, values are those used in the scans over the theoretical quantities after separation 
of experimental and theory errors. The first values of a and 7 are those uses in the comparison with CKMfitter and VTfi t . The 
second values, obtained from our global fit, are used in the studies comparing the effect of the inclusion of B(B + — > t + u). 



m t [GeV/c 2 ] 


m c [GeV/c 2 ] 


Am d [ps L ] 


Am s [ps *] 


173.3 ± 0.9 [13] 


1.275 ±0.025 [13] 


0.507 ± 0.004 [13] 


17.7 ±0.08 [13] 


Vcb 


v ub 


v ud 


K s 


(4.06 ±0.13) x 10" 2 [20] 
(4.09 ± 0.06 ± 0.11) x 10~ 2 [13] ^ 


(3.89 ±0.44) x 10" 3 [20] 
(4.15 ± 0.31 ± 0.48) x 10 -3 [13] a 


0.97425 ± 0.0002 [16] 


0.2208 ± 0.0039 [16] 




sin 2/3 


a 


7 


(2.228 ±0.0011) x 10" 3 [13] 


0.0682 ± 0.019 [20] 


(90.0 ±5)° [20] 
(84.0 ± 2.1)° 2, 


(76.0 ± 10)° [20] 
(78.8 ±4.4)° b 



a experimental and theory errors separated 
b Input from global fit (see below) 



TABLE II. QCD parameters used in the baseline fits with "statistical" and theory uncertainties, respectively. 



h. [MeV] 


fn a / fB d 


Bb 3 


B B JB Bd 


B K 


Ref. 


250 ±5.4 ±11 


1.215 ±0.012 ±0.015 


1.33 ±0.06 


1.05 ±0.06 


0.737 ± 0.006 ± 0.020 


[15] 


rice 


Vet 


Vtt 


rib 






1.39 ±0.35 


0.47 ± 0.04 


0.5765 ± 0.0065 


0.551 ±0.007 




f7] 



observables, we add terms in the x 2 , denoted generically 
by Mk, that include quark masses and B meson masses. 

We are careful to distinguish among different kinds 
of uncertainties. Observables with experimental er- 
rors only (statistical and systematic) are assumed to 
be Gaussian-distributed. Theoretical quantities such as 
lattice-derived QCD parameters, and inputs to V ub and 
Vcb, typically have two types of uncertainties. The first 
type of error is of a "statistical" nature, resulting from an 
input with statistical uncertainties or from Monte Carlo 
statistics in lattice calculations. We assume this error to 
be Gaussian-distributed and add corresponding terms to 
the x 2 i denoted by T n - The second type of uncertainty 
is a theory error with no known underlying statistical 
distribution. We therefore make no assumption as to 
the distribution of these errors, and instead perform a 
scan over a large range of plausible values, doing a \ 2 
minimization at each point. The scan includes the QCD 
parameters (f B , , f Bs / fB d , B Ba , B Bs / B Bd and B K ) and 
the CKM matrix elements V u b and V c b- The QCD cor- 
rections ry cc , r) ct and r\tt used in the determination of €k 
and rib that appear in the prediction of Amj also have 
theory errors. Although the fit methodology is able to 
scan over them, we do not do so, since this is unneces- 
sary at the current level of precision. We parameterize 
rice in terms of m c and a s @. Tables HI and [TXI summarize 
the input parameters for our baseline fits. We are prepar- 
ing a long article that includes a detailed discussion of all 



input values, provides a study of the correlations among 
the theory uncertainties and shows further results 0. 



III. COMPARISON TO CKMFITTER AND XJT fit 

To begin, we compare the performance of 
the scan method with that of CKMfitter and 
VJTfit using 19 input measurements (V u d, V us , 
V cb , V ub , e K ,Am d ,Am s ,sin2f3,a,~f, f Bs , B Bg , f B Jf Bd , 
B B jB Bd , B K , m t , m c , m Bd , m Bg ), choos- 
ing values specified for the book Physics 
of the B Factories [5j to fit 13 parameters 
(p, fj, A, A, f Bs , B Bs , f Bs /f Bd , B Bs /B bd ,B K ,m t ,m c , m Bd , 
m Ba ). In these fits we use the central values and mea- 
surement errors for a and 7 in the x 2 function. We 
compute mt(mt), which enters into the Inami-Lim 
functions [12| for Ato^, Ato s and €k, in the MS scheme 
from the pol e mass m t at three loop level for five 
quarks [9Hlll]. For fits having a P(x 2 ) probability 
> 32%, we plot la contours in the p — fj plane. For the 
central value, we select the fit with the highest P(\ 2 )- 
We take the ±lcr uncertainties from the maximum and 
minimum values of the envelope of all contours. We 
perform three different fits. In the first, we combine 
theory and experimental uncertainties, treating these 
as Gaussian. In the second, we scan over the theory 
uncertainties in f Bs ,B Bs ,BK- In the third, we separate 
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FIG. 1. Overlay of 68% CL contours in the p — r\ plane for 
the inputs of Physics of the B Factories with scanning over 
/b s , B Bs andB K • 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 



P 

FIG. 2. Overlay of 68% CL contours in the p — r\ plane for 
the inputs of Physics of B the Factories with scanning over 
/b, , B Ba , B K , V u b and V cb . 



TABLE III. Comparison of unitarity triangle parameters for different fitting techniques using inputs for Physics of the B 
Factories. The second and third columns show the results of fits by CKM fitter and UT/a. The fourth column shows the 
result of the fit if no scanning over theory parameters is performed. The fifth column shows the fit results if we scan over the 
QCD parameters /s s ,/s s //s d and Bk- The sixth column shows the fit result if we scan, in addition, over the theory errors 
associated with V u t and V c t- 



Parameter 


CKMfitter 


[5] 


Our fit 
no scan 


Scan i 
scan over 
/s s , }b s 1 /s d , B K 


nethod 

scan over 
/s s , }b s 1 'fB d , B K , V c b, V ub 


P 
V 

p n 

a [°] 
7[°] 


0.121 ±0.02 
0.349 ±0.012 
21.7 ± 1 
87.5 ±3.2 
70.9 ± 3.2 


0.125 ±0.022 
0.347 ±0.014 
21.6 ±0.8 
87.9 ±3.4 
70.4 ±3.4 


0.119 ±0.013 
0.353 ± 0.008 
21.9 ±0.3 
86.8 ± 2 
71.4 ±2 


1 16+ 034 
u.i±o_ .03i 

n qrq+0.020 

2i.8ig-» 

86.4j£i 

71.811;! 


1 ■H0+ ' 030 

r, oirr + O.OlS 
U.OUJ_q 

22.2±?: 4 
87.9l^: 3 3 
69.81^ 



theory uncertainties from experimental uncertainties 
in V u b and V c b and scan over the former. Figure [1] 
shows the contours in the p — fj plane for the second fit 
and Fig. [2] shows those for the third fit, i.e, the most 
complete scan. Table IIIII lists our results in comparison 
to those from CKMfitter and UT^. The three methods 
yield broadly similar results. Without scanning, the 
uncertainties in the unitarity triangle parameters are 
smaller than those from CKMfitter and UT^ t . Scanning 
only over the range of theory errors for the QCD 
parameters increases the uncertainties in the unitarity 
triangle parameters significantly. Scanning over V u b and 
V c b as well, the uncertainies are more than a factor of 
two larger than those from CKMfitter and UTfi t . The 
scan method finds acceptable fits in regions the other 



approaches do not explore, due to their assumptions as 
to the distribution of theory errors. This, of course, has 
an effect on the sensitivity of unitarity triangle fits to 
the contributions of physics beyond the SM. In view of 
the importance of searches for new physics in the flavor 
sector, it is worthwhile to use a procedure in making 
consistency tests of the unitarity triangle that treats the 
unknown statistical distribution of many theory errors 
with a minimum of assumptions. In this way claims of 
inconsistencies in unitarity triangle fits that violate the 
SM can be avoided. 

We also perform fits with the PDG12 V u b and V c b av- 
erages 13] that are listed in Table HI Since the V u b and 
V c b results from exclusive modes are significantly lower 
than those from inclusive modes, the PDG uses scaling 
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factors on the total errors of 2.6 and 2.0, respectively. In 
order to use the scan method in this case, we separate 
the experimental and theory uncertainties. In all subse- 
quent global fits of the CKM matrix, we scan over theory 
uncertainties in V u b, V c b, as well as those for fg a , JbJ fB d 
and Bk- Figure [2] shows the overlay of 68% confidence 
level contours in the p — fj plane of all accepted fits. 



v i 



E ft \ 


/ //''""' m V 

S /if —-f— 






: £ K / Y — 

iWvj I"" - ■ 

= I I ■ i I 1 I I 1 1 1 I I 1 ,*T 1 l-H 


lr i , ^ i , li 1 t ,1 , , i , , , tT>vj 


1 -0.8 -0.6 -0.4 -0.2 -0 0.2 0.4 0.6 0.8 

P 

3. Overlay of 95% CL contours in the p — fj plane i 



the baseline fit scan with 19 measurements without including 
B{B+ -> t+v). 
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FIG. 4. Overlay of 95% CL contours in the p — fj plane for 
the baseline fit scan with 19 measurements with inclusion of 
B(B+ -> r+v). 

IV. RESULTS OF THE GLOBAL FITS 
A. Fit results with direct inputs of a and 7 

To test the SM using the scan method, we relax the 
requirement for accepted fits to P(x 2 ) > 5%. If we have 
any fits that meet this criterion, the SM is deemed to 
be compatible with the data at the present level of the- 
oretical uncertainities. We form a contour correspond- 
ing to an accepted fit in any desired two-dimensional 
projection of parameter space. The contour consists of 
the boundary of the projection onto the plane of the 
multi-dimensional ellipsoid in parameter space satisfy- 
ing P(x 2 ) > 5%. We'll loosely refer to the contours as 
95% CL contours. Taking the extrema of such a contour 
along a parameter axis provides a nominal 95% CL inter- 
val for the parameter. The procedure actually undercov- 
ers somewhat, in the limit of no theoretical uncertainties. 



However, the effect on coverage is not more than a few 
percent for the cases encountered here. Figure [3] shows 
the overlay of 95% CL contours in the p — fj plane for 
all accepted baseline fits using 19 measurements to fit 13 
parameters, as in Section III. To obtain a range of val- 
ues for a given parameter, we take the extrema of the 
union of the contours attached to accepted fits. Table HVl 
shows the range of the unitarity triangle parameters thus 
obtained (hereinafter referred to as the 95% CL range). 

The B + — > t + v branching fraction [13], measured by 
BABAR [l8| and Belle [lj|, is very sensitive to contri- 
butions from a charged Higgs boson. The PDG av- 
erage B(B+ -> t+v) = (1.66 ± 0.33) x 10" 4 PHH 
is larger than the SM prediction of B(B+ — > t+v) = 
(1.2 ± 0.25) x 10~ 4 '21}. Figure H shows the results in 
the p — fj plane and Table IIVI summarizes the 95% CL 
ranges of unitarity parameters. Even with these high 
values of B(B+ — > t+v), we obtain a sizeable allowed 
p — fj region; there is no conflict with the SM. Belle has 
now presented a new measurement of B(B+ — > t+v) = 
(0.72ltg;!£(stat)±0.11 (sys)) x 10" 4 H3 that reduces the 
world average, and in turn reduces any perceived conflict 
in the SM with the other measurements of the unitarity 
triangle. 
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TABLE IV. The 95% CL ranges for unitarity triangle parameters from our baseline fit scans without and with the inclusion of 
B(B + — > t + u). The values of a and 7 are, in these cases, computed from the values of p and fj. 



Parameter 


P 


n 


p n 


a n 


7[°] 


Scan without B + — >• t + u 
Scan with B + — > t + v 


0.058 - 0.181 
0.085 - 0.159 


0.324 - 0.395 
0.334 - 0.377 


20.8 - 23.7 
21.2 - 22.9 


77.1 - 95.9 
76.9 - 92.9 


62.3 - 81.0 
65.2 - 76.9 
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FIG. 6. Overlay of 95% CL contours in the a — j3 plane for 
B ->■ PP, B — >• PV and B -> VV modes. 



B. Fit Results using Individual Measurements of 
Branching Fractions and CP Asymmetries of Various 
Decays 

We also perform fits in which, instead of treating the 
values of a and 7 as inputs, we directly include the mea- 
surements that determine them. This allows us to de- 
termine the correlations in the extraction of the var- 
ious unitarity triangle angles from the fit. We omit 
B(B + — > t + v) in these fits; this has little effect on the 
results. 

For the a determination, we include all measured 
branching fractions and CP asymmetries in B — > PP, 
B ->• PV and B — > VV modes in the fit. Follow- 
ing the Gronau-Rosner approach (23| . we parametrize 
amplitudes in terms of tree, color-suppressed tree, pen- 
guin, singlet penguin, VF-annihilation/VF-exchange, elec- 
troweak and color-suppressed electroweak diagrams (up 




1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

FIG. 7. Overlay of 95% CL contours in the a — ft plane 
from fits including the B -> PP, B -¥ PV, B -¥ VV and 
B — > Pa\ modes. The vertical band shows the 68% CL range 
for P obtained from sin 2/3 measurements [2[| . 

to A 3 beyond leading order). Thus for order A 3 , we con- 
sider SU(3) corrections for the tree, color-suppressed tree 
and penguin diagrams. Since calculations of branching 
fractions use the B + and B^ life times, we add \ 2 terms 
for them in the fit. 

For the 7 determination, we include measured branch- 
ing fractions and CP asymmetries for B + — > D^K + 
and B + — > DK* + decays analyzed in the Giri- 
Grossman-Soffer-Zupan (GGSZ) j24|, Gronau-London- 
Wyler (GLW) [H and Atwood-Dunietz-Soni (ADS) 
methods (26|. We also include branching fractions and 
CP asymmetries for B + D^n + decays analyzed in the 
ADS method. Here, the observables are calculated in 
terms of b — > cus(d) and b — > ucs(d) amplitudes. We 
also include time-dependent CP asymmetries in B° — >• 
£)(*)+7r- an d B° -4 DW+ p~ decays that determine 
sin(2/3 + 7). Thus, the total number of measurements 
in the global fits increases to 230; the fit has 104 parame- 
ters. Figure [S] shows the 95% CL contours of all accepted 
fits in the p — fj plane. The allowed region is smaller than 
that for the baseline fits. Nonetheless, this result shows 
that the 230 measurements are in good agreement with 
the SM. This procedure accounts for possible correlations 
between a, 7 and the other Wolfenstein parameters (e.g., 
P). 



V. DETERMINATION OF a AND 7 

In addition to the global fits of the CKM matrix, we 
perform separate fits that determine the unitarity trian- 
gles a and 7. 



6 



A. Determination of a 

For the a determination we combine all measured 
branching fractions and CP asymmetries in B — > PP, 
B — > PV and B — > VV modes. We first perform sep- 
arate fits for each class of decays. As in Section IV we 
parameterize the observables in terms of amplitudes, fol- 
lowing the Gronau-Rosner method. We include Jb s and 
fB s /fB4 m the fit j but we do not scan over them, since 
these parameters appear only in the IV-exchange and 
IV-annihilation diagrams that are at order A 3 ; any vari- 
ation in these parameters is absorbed by adjusting the 
magnitude of the IV-annihilation/IV-exchange diagrams 
leaving a and the remaining parameters unchanged. We 
perform a single fit for each class of decays and plot 
95% CL a — (3 contours. Figure [5] shows the a — (3 
contours at 95% CL for the fits to B -> PP, B -> PV 
and B — > VV separately. The contours in the a — (3 
plane, which are rather wide, encompass the world aver- 
age j3 = (21.4 ± 0.8)° measured using b — > ccs modes. 

We next perform a combined fit of all measurements 
in B -> PP,B -> PV and B — > VV and B ->■ Pa x 
modes to extract a. We use 181 measurements [20J to 
determine 94 parameters. Figure [7] shows the 95% CL 
contours in the a — (3 plane, which delineate a region 
that is considerably smaller than the contours of the fits 
for the individual decay classes. This is due to the fact 
that the B -t PP, VV (B -> PV) modes yield larger 
(smaller) central values of j3 than the sin 2/3 measurement 
in b ccs decays (see Table fV|) . The fit probability is 
P(X 2 ) = 9.6%. Table ED lists all fit results for a and 
(3. These results show that all measurements in B — > 
PP, B ^ PV and B -> VV and B -> Pa x modes are 
consistent with the SM description and no new physics 
amplitudes are required. 



B. Determination of 7 

For the 7 determination, we use branching fraction and 
CP asymmetries of B + -> D^K + and B+DK*+ de- 
cays analyzed in the GGSZ H, GLW g| and ADS [26| 
methods. We also include branching fractions and CP 
asymmetries of B + — >• D^n + decays analyzed in the 
ADS method and time-dependent CP asymmetries in 
B° -> dM+tt- and B° -> flW+ p~ decays that de- 
termine sin(2/3 + 7). We separate the CKM factors, 
\V us V: b \/\V cs V c * b \ and \V ud V: b \ / \V cd V c * b \ from the ratio of 
b — > u to b — >• c amplitudes. We include the ratios 
\V US /V ud \ and \V ub /V cb \ in the fit, scanning over \V ub \ and 
\V cb \. Since the predictions contain a product of CKM 
factors and ratios of amplitudes, it is necessary to con- 
strain the CKM factor in the fit to obtain sensible values 
for the amplitude ratios and CKM factors. We use 56 
measurements [20] to extract 19 fit parameters, scanning 
over the constraint V ub | / 1 V cb | . The probabilities for these 
fits range between P(x 2 ) = 10.4% and P{x 2 ) = 14.4%. 
Figure [5] shows the resulting contours at 95% CL in the 



7 — (3 plane. Table [V] lists the fit results. Again, these 
56 modes are well-described within the SM and no new 
physics amplitudes are required. 




FIG. 8. Overlay of 95% CL contours in the 7 — ft plane for 
fits including the B -¥ DK M , DK* , D W 7r, Dp modes. The 
vertical band shows the 68% CL j3 region obtained from sin 2/3 
measurements 12011 . 



VI. CONCLUSION 

The three fitting approaches: CKMfittcr, UT^ and 
the scan method yield similar central values for p and fj 
when presented with identical inputs. If we employ the 
scan method procedure of scanning over theory uncer- 
tainties in the QCD parameters, and V ub and V cb as well, 
making no assumption as to the distribution of theory 
errors, however, the allowed region in the p — fj plane is 
significantly larger. As a result, therefore, in the scan 
method measurements of the unitarity triangle provide 
good agreement with the SM; there is no need at this 
juncture to invoke new physics. The primary reason for 
the enlarged size of the allowed region in the p — f\ plane 
is the large theory uncertainties in V cb , V ub and the QCD 
parameters, coupled with the lack of knowledge about 
their distribution. As there is no principled frequen- 
tist method of ascertaining the distribution of theory er- 
rors, absent a Bayesian belief in the relevance of a Gaus- 
sian prior, the errors in p and fj resulting from the scan 
method determination are obtained with a minimum of 
assumptions. They are thus less likely to result in spuri- 
ous indications of a breakdown of the SM. 

For this reason, using the scan method, we find no 
tension with the SM even when we include the current 
PDG value of B(B+ -> t+u) = (1.66 ± 0.33) x 10~ 4 ; the 
scan yields global fits consistent with the SM at 95% CL. 
When we include the recent Belle result as well, the al- 
lowed region becomes even larger. Using all measured 
branching fractions and CP asymmetries of B — >• PP, 
B -> PV, B -4 VV modes and B+ -4 D^K+,B -> 
DK* + modes that are sensitive to a and 7, respectively, 
we see a reduced allowed region in the p — fj plane. From 
separate fits of branching fractions and CP asymmetries 
in these modes, we determine a — (3 and 7 — (3 contours. 
Though a measurements agree with each other and (3 re- 
sults are consistent with sin 2/3 from b — > ccs modes, the 
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TABLE V. Measurements of a. ft and 7 from fits of branching fractions and CP asymmetries in B — > PP, B — > PV, B — > VV 
decays and in a combination of all modes combined plus B — > Pai decays (a) and B —>■ D^K{n) + B — > DK*(p) decays (7). 





B — > PP 


B^ PV 


B VV 


B modes combined 


B -» D^K{ir) + B -S- DK*(p) 


7[°] 
P [°] 


86.5 ±3.4 

25.8±|;2 


80.8 ±4.0 

i8.3l|:5 


81.4 ±5.2 
27.5l^; 4 7 


84l 2 2 ;J 

20.8li;^ 


78.811* 
22.81 2.1 



correlations among Wolfenstein parameters in the differ- 
ent measurements appear to be important. The values 
of a, determined from a fit to all measured branching 
fractions and CP asymmetries of B — > PP, B — > PV, 



B — > VV and B — > Pa± modes, and 7, extracted from 
a fit to B+ -> D^K+,B -> DK*+, B -> £>«7T and 
B° — > D + p~ modes, agree with the SM expectations. 
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